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A two-dimensional semi-empirical model of pulsed inductive thrust efficiency is devel- 
oped to predict the effect of such a geometry on thrust efficiency. The model includes 
electromagnetic and gas-dynamic forces but excludes energy conversion from radial mo- 
tion to axial motion, with the intention of characterizing thrust efficiency loss mechanisms 
that result from a conical versus a flat inductive coil geometry. The range of conical pulsed 
inductive thruster geometries to which this model can be applied is explored with the 
use of finite element analysis. A semi-empirical relation for inductance as a function of 
current sheet radial and axial position is the limiting feature of the model, restricting 
the applicability as a function of half cone angle to a range from ten degrees to about 60 
degrees. The model is nondimensionalized, yielding a set of dimensionless performance 
scaling parameters. Results of the model indicate that radial current sheet motion changes 
the axial dynamic impedance parameter at which thrust efficiency is maximized. This shift 
indicates that when radial current sheet motion is permitted in the model longer character- 
istic circuit timescales are more efficient, which can be attributed to a lower current sheet 
axial velocity as the plasma more rapidly decouples from the coil through radial motion. 
Thrust efficiency is shown to increase monotonically for decreasing values of the radial dy- 
namic impedance parameter. This trend indicates that to maximize the radial decoupling 
timescale should be long compared to the characteristic circuit timescale. 


Nomenclature 


c 

capacitance (F) 

Rpi Re 

plasma, external resistance (fi) 

E 

energy (J) 

• coil i ' coil 

coil minor, average radius (m) 

F 

force (N) 

Ti 

upstream temperature (K) 

h,h 

coil, plasma current (A) 

u 

shock velocity (m/s) 

k 

Boltzmann’s constant (J K -1 ) 

V,Fo 

voltage, initial capacitor voltage (V) 

M)i Lq, L tot 

initial, coil, total inductance (H) 

v r ,v z 

radial, axial slug velocity (m s -1 ) 

L* 

nondimensional inductance 

z 

axial location (m) 

Icoil 

coil length (m) 

Zo 

axial decoupling distance (m) 

M 

mutual inductance (H) 

a 

axial dynamic impedance parameter 
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M 

Mach number 

7 

specific heat ratio 

77li , TTlbit 

ion mass, slug mass (kg) 

<t> 

radial dynamic impedance parameter 

Pl,P2 

upstream, downstream pressure (Pa) 


critical resistance ratios 

r, f 

radial, average radial location (m) 

E 

gas-dynamic pressure parameter 


I. Introduction 

I T is desirable to extend the lifetime and increase the reliability of an in-space propulsion system as much 
as possible since maintenance or replacement of that system becomes particularly challenging once it 
has been placed into orbit. Reducing the size and mass of the propulsion system, including the propellant 
required to complete a mission, can permit an increase in the amount of payload as a percentage of total 
vehicle mass. The high values of specific impulse associated with electric propulsion (EP) reduces the amount 
of propellant needed for a given mission relative to other conventional propulsion systems. 

Pulsed inductive plasma thrusters [1-3] are spacecraft propulsion devices in which electrical energy is 
capacitively stored and then discharged through an inductive coil. The thruster is electrodeless, with a time- 
varying current in the coil interacting with a plasma covering the face of the coil to induce a plasma current. 
Propellant is accelerated and expelled at a high exhaust velocity (O (10 — 100 km/s)) by the Lorentz body 
force arising from the interaction of the magnetic field and the induced plasma current. 

Thrusters of this type possess many demonstrated and potential benefits that make them worthy of 
continued investigation. The electrodeless nature of these thrusters eliminates the lifetime and contamination 
issues associated with electrode erosion in conventional electric thrusters. Also, a wider variety of propellants 
are available for use when compatibility with metallic electrodes is no longer an issue. Pulsed inductive 
accelerators have been successfully operated using propellants like ammonia, hydrazine, and CO 2 , and there 
is no fundamental reason why they would not operate on other propellants like H 2 O. It is well-known that 
pulsed accelerators can maintain constant specific impulse, J sp , and thrust efficiency, r/t , over a wide range 
of input power levels by adjusting the pulse rate to maintain a constant discharge energy per unit pulse. It 
has also been demonstrated that a pulsed inductive thruster operating in or near the regime of optimum 
dynamic impedance matching can operate at a relatively constant thrust efficiency over a wide / sp range. 
Thrusters in this class have operated at high energy per pulse, and by increasing the pulse rate they offer 
the potential to process very high levels of power using a single thruster. 

It has been found [1,2,4] that propellant utilization inefficiency can be a disadvantage for pulsed inductive 
thrusters with a flat inductive coil geometry. One proposed solution is the use of a nozzle downstream of the 
accelerating coil that injects propellant upstream onto the face of the inductive coil. This solution, however, 
presents the additional disadvantage of placing a physical body in the exhaust path of the propellant, 
adversely altering the propellant trajectory as the exhaust impacts and, over time, erodes the propellant 
injection structure. 

A possible alternative solution to this problem is to alter the inductive coil geometry to better confine 
the neutral propellant within the region of interest where the processes of current sheet formation and 
acceleration occur. One alternative to a flat coil is a conic section or conical frustum. A flat, disc-shaped 
inductive coil can be considered as a conical inductive coil with a half cone angle of 90°. As the half cone 
angle is decreased the fraction of axial electromagnetic force on the plasma current sheet decreases as a radial 
electromagnetic (pinching) force appears and increases. The ability to predict the effect of coil geometry 
on propellant utilization efficiency and plasma current sheet motion would permit a trade-off study between 
cone angle and efficiency and aid in the physical understanding and design of this type of thruster. 

A model [5] relating two-dimensional plasma current sheet motion and inductive coil geometry has been 
developed based upon a two-dimensional semi-empirical inductance relation incorporated into a well-known 
acceleration model [2]. The semi-empirical model yielding inductance as a function of plasma position has 
been compared to experimental data obtained for a limited range of coil geometries and simulations spanning 
a broader expanse of the parameter space. In this paper we present the semi-empirical model of inductance 
as a function of coil angle and current sheet position and show finite element results that support the model’s 
validity. The entire two-dimensional acceleration model is then nondimensionalized and the relevant scaling 
parameters are discussed in terms of their effect on thrust efficiency. 
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II. Two-Dimensional Inductive Accelerator Model 


Pulsed inductive thrusters are commonly modeled as a circuit [2], shown in Fig. II, where /i is the current 
flowing in the driving circuit, I 2 is the current flowing in the plasma current sheet, C is the capacitance of 
the capacitor, M is the time- varying mutual inductance between the driving coil and the current sheet, Lq 
is the initial (parasitic) inductance, Lq is the accessible coil inductance, R p is the resistance of the plasma, 
and R e is the resistance in the driving circuit. 


R e Lo Switch 



Figure 1. Lumped element circuit model of a pulsed inductive thruster inductively coupled to a plasma (taken 
from Ref. [2]). 


R e Lo Switch L c — M L c — M 



Figure 2. Equivalent circuit of a pulsed inductive thruster inductively coupled to a plasma used to derive a 
set of governing equations to model thruster performance. 


An equivalent circuit, shown in Fig. 2, can be drawn and a set of equations can be written to model 
thruster performance in terms of the electrical parameters. For this circuit, where V is the voltage on 
the capacitor as a function of time (initially at Vq), the following equations follow from the application of 
Kirchoff’s law: 


dli LcV — LcReh ~ MRpI -2 + (L c h + MIi) 
dt Lc {L 0 + Lc) — M 2 


(1) 


d/2 

dt 


M 


dli 

dt 


dM 

dt 


II — Rpl2 


Lc 


dV _ _/i 
dt C 


( 2 ) 

( 3 ) 


Neglecting the resistive elemenets and adding the series and parallel lumped inductive elements shown 
in Fig. 2 gives the total inductance, 


Ltot — Lq + Lc 


M 2 

Lc 


( 4 ) 


It can be seen from this equation that the total inductance changes as a function of time due to the mutual 
inductance term, which varies due to the movement of the current sheet. An expression for the mutual 
inductance as a function of current sheet axial and radial position has been empirically determined as: 


L to t(r,z) 


Lo + Lc 


exp {—z/zq) 



( 5 ) 
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where f is the average radial position of the current sheet, r co u is the average radial position of the coil, z is 
the axial displacement of the current sheet, N is a fit parameter found to be close to 2 for all cases studied 
here, and Zq is the axial decoupling distance. This expression is set equal to the previous expression for total 
inductance (Eq. 4) and solved for the mutual inductance as a function of the axial and radial separation 
distance between the driving coil and the current sheet: 

M = L c exp(^) ( 6 ) 

\ 2-2^0/ \T coil / 

of which the time derivative is, 


dM 

dt 


Lc N _ k _ i 

r 2 1 

N o 

T coil 2 


dr 

dt 


exp(— 2 /zq) 


L C 

2zq 


exp (-z/2zq) 


dz 

dt 



(7) 


The equation of motion for the current sheet can be written using Newton’s second law with the thrust- 
producing force arising from the magnetic pressure between the driving coil and the current sheet. This 
force in the arbitrary i "'-direction is given as: 


Fi = 


I 2 dL 
2 dxi ' 


(8) 


This force is opposed in the radial direction by a gas-dynamic pressure force P 2 that increases as the current 
sheet moves toward the thrust axis and is calculated assuming the current sheet acts as a normal shock wave: 


^ = 1 +^ 
Pi 7+1 


[M 2 - 1] 


(9) 


where P 2 is the pressure of the gas downstream of the shock, Pi is the pressure upstream of the shock, 7 is 
the ratio of specific heats (taken here to be 5/3), and A4 is the local Mach number upstream of the shock: 


m = u,m ( 10 ) 

V TOj 

where u is the shock velocity in the given direction, k is Boltzmann’s constant, T± is the temperature of the 
gas upstream of the shock, and m t is the mass of an ion. 

Using the functional dependence of inductance on both axial and radial position, scalar momentum 
equations can be written for both axial and radial motion: 


dv z 

dt 


LcA 

2 z 0 


exp 




( 11 ) 


dv r 

dt 


P 2 27t? 7 co ,: J - 


L c llN 


2 r* 


-N 


exp 


coil 


— M 

Zo 


N - 1 


/mbit. 


(12) 


where mbit, is the total propellant mass in the current sheet, v z is the axial current sheet velocity, and v r 
radial velocity. The axial 2 and radial r position are related to the axial and radial velocities as: 


is 


dz 

dt 


= v z 


(13) 



The time-evolution of the pressure P 2 is given by the time-derivative of Eq. 9, 


(14) 


dP 2 Pi 2 y mi dv r 

(15) 

bringing the total number of first-order coupled ODEs in the system to nine. 

Of the nine equations (Eqs. 1, 2, 3, 7, 11, 12, 13, 14, and 15), six rely on the semi-empirical relation 
for coil inductance as a function of current sheet location. It should be noted that no attempt is made here 
to model energy conversion from radial current sheet motion to axial momentum. In addition, this study 
focuses only on thrust efficiency, making no attempt to model factors such as propellant utilization efficiency 
that could be a function of coil geometry and would impact the total thruster efficiency. 
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III. Finite Element Analysis 


The inductive coupling between coils and current sheets of various geometries were simulated using finite 
element analysis to explore the geometric range of validity of Eq. 5. First, the coil inductance Lc and 
decoupling length Zq were calculated by curve fitting Eq. 5 to the simulated inductance as a function of 
axial current sheet position with a constant current sheet radial position roughly equal to (r co u) . Then, the 
inductance was calculated for various current sheet radial and axial locations, the results of which are shown 
in Figs. 4, 5, and 6 along with inductance calculations using Eq. 5 where the simulation data are plotted as 
markers and Eq. 5 is plotted as a red line. 



Figure 3. Diagram of the inductive coil geometry. 


The coil geometry is shown and labeled in Fig. 3. The geometry is designated in the lower left hand 
corner of Figs. 4-6 where the number is the half cone angle and a suffix of “S” , “M” or “L” refers respectively 
to a short, 5 cm, medium 6-8 cm, or a long, 10 cm coil length l co u. All coils had a minor radius r co u of 4 cm, 
and the only coil with l co u = 0 had a half cone angle of 90°. Agreement between the finite element results 
and Eq. 5 is good for half cone angles between 20°-55°, as shown in Fig. 5, but begins to breakdown outside 
this range, as shown in Figs. 4 and 6. 





Figure 4. Finite element results (represented by markers) and the semi-empirical model (represented by a red 
line) for the inductance as a function of average current sheet radial position. Coil geometry is designated in 
the lower-left corner of the plot. 


For coils with a half cone angle below 12°, it is not the radial part of the model that begins to fail but 
rather the axial part. The inductance as a function of axial current sheet location is fit more accurately by 
an error function than an exponential, as shown in Fig. 7. 

Plots of the inductance as calculated using the error function in place of the exponential in Eq. 5 are 
compared in Fig. 8 to inductance calculations using Eq. 5. Below half cone angles of 12°, the discrepancies 
in the results of two semi-empirical equations for inductance increase with decreasing half cone angle, with 
the exponential function providing a poorer fit to the finite element inductance results. 
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Figure 5. Finite element results (represented by markers) and Eq. 5 (represented by a red line) for the 
inductance as a function of average current sheet radial position. Coil geometry is designated in the lower-left 
corner of the plot. 





Figure 6. Finite element results (represented by markers) and Eq. 5 (represented by a red line) for the 
inductance as a function of average current sheet radial position. Coil geometry is designated in the lower-left 
corner of the plot. 


IV. Nondimensional Solution Approach 


The following substitutions can be made to nondimensionalize the equation set from Section II following 
Refs. [6,7]: 


It = 


1 Lc 


Vo V c 


II 


V 

v* = — 

Ro 

* VLoC 
V z = ~—Vz 

z 0 

* VL^C 

v r = V r 

T coil 

t 

s/UC 


t* = 


T* — 

lr> 


1 Lc 


Vo V C 


M* = ^~ 
Lc 

* z 
z* = — 
Zo 


r = 


(16) 


p* 


T coil 

P 


Pi 


Applying these substitutions to Eqns. 1, 2, 3, 7, 11, 12, 13, 14, and 15 yields the following nondimensional 
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Figure 7. Comparison between the functional form of fits for the inductance as a function of current sheet 
axial location. Exponential fits are shown in red and error function fits are shown in blue. 


equation set: 


dl* 

dt* 

dll 

dt* 

dV* 

Ft* 

dM* 

dt* 

dv* 

Ft* 

dr* 

dt* 

<K 

dt* 

dP* 

Ft* 


7 7 \/f* 

L*V* + ( M*I{ + I*) - i\> 2 L*I*M* - il)iL*I{ 


/( 


L* + 1 - M* 


M* — + 7? InL 1p2 

dt* 1 dt* 2 r 


-I* 

N * . z* 1 

—r v r exp ( — j —v v z exp(- — ) 

XP*r* - (j>lf r* N ~' exp(-;z*) 


V 


* 

r 


a I* 2 r* N exp(— z*) 


dv* 

Ft* 


( 17 ) 


The terms a , ipi, if) 2 , 4>, A, and S are the relevant nondimensional scaling parameters of the system, and 
are defined as: 


A = 


V 0 2 C 2 L c 
2m bit zl 
C_ 
Lq 

LqC Pi2jTtI co h 


i >2 = R P \ -r~ 


L* = 


2f Tlbit 

Lo 


Lc 


F = F 


^ _ VqC 2 L c 

‘^TTlbit'y'coil 
4y rrii 1 
7+1 7 fcTi r 2 coil L 0 C 


( 18 ) 
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Figure 8. Comparison of the accuracy of two functional fits for the inductance as a function of current sheet 
axial and radial locations. Exponential fits are shown in red and error function fits are shown in blue. 


V. Model Results and Discussion 

The parameters if) \ and if ) 2 are critical resistance ratios, affecting the damping of the current in the driving 
circuit, and L* is a ratio of the initial (or inaccessible) inductance to the accessible inductance, representing 
a measure of the fraction of input electromagnetic energy that can be transferred to acceleration of the 
current sheet. The physical interpretations of these three parameters remains unchanged as compared to 
those of previous studies of pulsed inductive plasma thrusters with flat inductive coil geometries, and more 
information on them can be found in Refs. [6,7]. 

A appears as a new parameter in this model, and affects the force opposing radial current sheet motion 
caused by the gas-dynamic pressure acting in the positive radial direction along the inner surface of the 
current sheet. 3 is another new parameter, and is a measure of the growth rate of this gas-dynamic pressure 
as the current sheet shock front undergoes radial acceleration. For the following calculations, ip\ = 0.05, 
V > 2 = 0.01, L* = 0.18, A = 7.0 • 10" 7 , and 3 = 2.3 • 10 4 . 

The parameter a can be expanded into a product of physically meaningful ratios [ 6 ] : 

C 2 V 2 L c = 1 CV 2 / 2 ^ / WZg7 \ 2 
2 m bit zl 87 r 2 m bit v 2 /2 \ L 0 /L z ) 

where L z = v z L' z and L' z = Lq/zq. The second ratio on the right hand side represents the inverse of thrust 
efficiency, the third term represents a ratio of the initial inductance in the pulse circuit to the coil inductance, 
which places a limit on the achievable thrust efficiency through the Lovberg criterion [ 8 ] . The fourth ratio 
represents a balance between the characteristic ringing time of the driving circuit and the characteristic time 
over which the current sheet remains electromagnetically coupled to the driving circuit. In other words, the 
parameter a is a measure of the dynamic impedance match between the driving circuit and the axially-moving 
current sheet. It has been shown that there exists an optimum a for which r) t is maximized, corresponding 
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to a pulsed inductive plasma thruster that is dynamically impedance matched [9] . 
Thrust efficiency can be related to a by the following relation: 


For systems constrained to axial motion (zero radial motion) the optimum value of a for maximum r/t is 
between one and three [10]. When the current sheet is axially accelerated at a fixed radial position that is less 
than r co u (v* = 0 and r* = AW t* where 0 < A < 1), the absolute value of r/t decreases and the maximum 
r/t occurs at a higher value of a, as shown in Fig. 9. As shown in the figure, the maximum achievable r/t is 
lower as the sheet is displaced closer to the centerline. 



Figure 9. Calculated thrust efficiency as a function of a for different radial current sheet displacements. 

A possible explanation for the shift in optimum values of a is that as the current sheet displaces radially, 
it slows the axial acceleration of the current sheet. This causes a longer current sheet residence time in 
practice, but that information is not contained in a. Consequently, a longer characteristic circuit time, 
which will increase the value of a, represents a better dynamic impedance match. 

The parameter <f) can be similarly expanded: 

C 2 Vq L c __ 1 C% 2 / 2 / 2t tVL^C \ 2 

2 m. blt rl oil 8n 2 m b itV%/2 \ L 0 /L r ) 

where L r = v r L' r and L' r = Lc/r co u. The second term does not represent thrust efficiency, as in the 
expansion for a, but rather the inverse of the fraction of total energy that is converted into radial kinetic 
energy. The fourth term represents a balance between the characteristic circuit time and the characteristic 
time over which radial (not axial) current sheet motion causes electromagnetic decoupling from the driving 
circuit. 

While an efficient coupling of energy into the axial direction is beneficial to thrust efficiency, efficient 
coupling of energy into the radial direction has a generally detrimental effect on thrust efficiency unless 
some means exists to efficiently convert the radial motion to axial thrust. If the radial dynamic impedance 
mismatch were caused by a small timescale of radial current sheet decoupling compared to the timescale of 
the driving circuit, the current sheet would decouple before the driving circuit had an opportunity to transfer 
initial circuit energy into axial current sheet motion. Therefore, thrust efficiency should be optimized for a 
condition where the timescale over which the current sheet radially decouples from the driving circuit is long 
compared to the characteristic time of the driving circuit, leading to faster decoupling due to axial motion 
instead of radial motion. This would indicate an increase in r] t for decreasing values of </>, a trend that can 
be seen in the plot of rjt versus <j> for a = 0.6, shown in Fig. 10. 

A contour plot of the combined effects of the two dynamic impedance parameters is shown in Fig. 11. As 
</> increases, the thrust efficiency decreases for all values of a due to faster current sheet decoupling from the 
driving circuit, resulting in lower achievable axial propellant velocity. The efficiency is optimized at higher 
values of a as 0 is increased, in agreement with Fig. 9. 
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Figure 10. Model results of thrust efficiency as a function of <f> for the unconstrained case with a = 0.6. 
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Figure 11. Contour plot of thrust efficiency as a function of a and <f>. 

VI. Conclusions 

In this paper we presented a two-dimensional semi-empirical acceleration model for pulsed inductive 
plasma thrusters. The semi-empirical aspects of the model were validated through a combination of mea- 
surements and finite element modeling of the inductance of coil geometries with a minor radius of 4 cm 
where the half cone angles were between 20°-55° and the coil lengths ranged from 5-10 cm. The model was 
nondimensionalized and two of the nondimensional scaling parameters, namely the axial and radial dynamic 
impedance parameters, were discussed. Contour plots of thrust efficiency as a function of these two param- 
eters were presented. Trends in these plots indicate that as 4> is increased, thrust efficiency is optimized at 
greater values of a. This is explained by a slower axial current sheet acceleration and consequently a longer 
axial acceleration timescale due to radial current sheet displacement, which results in a better dynamic 
impedance match occurring for a longer characteristic circuit time. In the present model with no mechanism 
to convert radial kinetic energy to directed axial thrust, efficiency is maximized for smaller values of the 
radial dynamic impedance parameter, indicating that higher axial current sheet velocity is achieved when 
the current sheet decouples slowly in the radial direction compared to the characteristic time of the driving 
circuit. 
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